Coupling-Matrix Approach to the Chern Number Calculation in Disordered Systems 
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The Chern number is often used to distinguish between different topological phases of matter 
in two-dimensional electron systems. A fast and efficient coupling-matrix method is designed to 
calculate the Chern number in finite crystalline and disordered systems. To show its effectiveness, 
we apply the approach to the Haldane model and the lattice Hofstadter model, the quantized 
Chern numbers being correctly obtained. The disorder-induced topological phase transition is well 
reproduced, when the disorder strength is increased beyond the critical value. We expect the method 
to be widely applicable to the study of topological quantum numbers. 
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I. INTRODUCTION 

In the past thirty years, the condensed matter physics 
community has been fascinated by topological phases of 
matter, for instance, the integer quantum Hall effect, 1 
the fractional quantum Hall effect^ the quantum anoma- 
lous Hall effect,— the quantum spin Hall effect,— £ and the 
three-dimensional topological insulatorsi&l These topo- 
logical states of quantum matter are usually distin- 
guished by use of some global topological quantum num- 
bers^ rather than certain local order parameters. The 
topological aspect of the integer quantum Hall effect 
with periodic potentials was first discussed by Thouless, 
Kohmoto, Nightingale, and Nijs (TKNN)^ In their fa- 
mous work, a topological expression for the Hall conduc- 
tivity was given by the Chern number over the magnetic 
Brillouin zone. Their result was then generalized to the 
fractional quantum Hall effect^ For the quantum spin 
Hall systems, with the extension of the idea, the well- 
defined spin Chern number can be used to characterize 
trivial and non-trivial bulk band topology 11 ' 12 

While simplification exists for pure systems^, calcula- 
tion of the Chern number in the presence of disorder is 
usually based upon the integral of partial derivatives of 
electron wave functions over the boundary phasesj 10 ' 14 : 15 
Numerical implementation involves hundreds of times 
of exact diagonalization for a given disorder configura- 
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which is very time-consuming even for nonin- 
teracting electron systems. Recently, several different ap- 
proaches for the Chern number computation in real space 
have been suggested in the literature. Kitaev^ proposed 
a real space Chern number formula for a lattice model 
in terms of traces of the coordinate operator and projec- 
tion operator. A recent work similar to Kitacv's one was 
numerically realized by Bianco and Resta 19 , and they 
generalized the result to the systems with open bound- 
ary conditions. In the presence of disorder, however, 
the idea of supercells,— namely, a periodic duplication 
of the actual system is needed, which greatly increases 
the computation time. Based upon the noncommutative 
Chern number theory, Prodan et alM- proposed an effi- 
cient method to calculate the Chern number of disordered 



systems, in which the procedure of the exact diagonaliza- 
tion is greatly simplified. However, it involves somewhat 
complicated multiple commutators between the coordi- 
nate operators and projection operators. In some cases, 
the Chern number can be extracted indirectly from the 
transport coefficients, which can be relatively easily cal- 
culated from the Kubo formula^ or the scattering ma- 
trix.— The latter approach 22 is based upon the bulk-edge 
equivalence, and is also applicable to the study of the Z2 
topological insulators. In some other cases, direct calcu- 
lation of the topological number of bulk wavefunctions 
is often needed and sometimes irreplaceable. Therefore, 
development of efficient numerical approaches to direct 
calculation of the Chern number is highly desirable. 



In this work, we provide an alternative way to calcu- 
late the Chern number, in which only one time exact 
diagonalization for the actual system is needed without 
loss of accuracy. A transparent coupling-matrix formu- 
lation will be given, from which the Chern number can 
be very efficiently computed, compared with the existing 
approaches. To show its effectiveness, this approach is 
applied to both the Haldane model and the Hofstadter 
model. The calculated Chern number is found well quan- 
tized provided the Fermi level lies within the energy gap, 
even when the sample size is not very large. The topo- 
logical phase transition from the quantum Hall insulator 
to an ordinary insulator can be determined based upon 
the calculated Chern number, and the obtained critical 
disorder strength is in good agreement with the result 
previously obtained from the Hall conductivity calcula- 
tion. 



In the next section, we present the new approach of cal- 
culating the Chern number. In Sec. Ill, for the Haldane 
model, it is shown that the present approach works well 
for both crystalline and disordered systems. In Sec. IV, 
we apply the approach to the lattice Hofstadter model 
and the calculated results are in accordance with already 
existing results. The final section is a summary. 
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II. METHOD 

We now consider a two dimensional (2D) lattice with 
N = L x x L y unit cells. We use r = (a;, y) with x and y 
as integers to index the position of a unit cell. We can 
define twisted boundary conditions for a single-particle 
wavefunction of the system <pg(x + L x ,y) — e l6 *Lpg(x, y) 
and <pg(x,y + L y ) — e l6y Lpg(x,y) with 6 = (8 x ,8 y ) and 
< 8 x ,8 y < 27r. The wavefunction is in general a row 
vector in the space of inner degrees of freedom, such 
as spin and sublattice. The system is assumed to have 
M electrons, and in the ground state the wavefunctions 
of the M occupied single-particle electron states are de- 
noted by <p™(r) with m — 0, 1, • • • M — 1. The many-body 
wavefunction of the ground state ^({rj}) is the Slater 
determinant of the single-particle wavefunctions (p™^), 
where r< with i — 0, 1, • • • M — 1 is the coordinate of the 
j-th electron. The Chern number of the gound state is 
given byi&Iiiis 

C=^-. I dO{V e ^e\ x |V e * e ) , (1) 

with 8 — (8 x ,8 y ) and Tg denoting the (8 x ,8 y ) space, 
which is essentially a torus. 

We transform the single-particle eigenstate from real 
space to momentum space through a Fourier transform 
(FT) 

^(r l ) = ^^F™(k i )e' ;k '^ . (2) 

The twisted boundary conditions require that the mo- 
menta take only the discrete values kj = k^ -* + q, where 
k[ 0) = f£) with < n < L x and < I < L y , and 

q = (j^ij^-)- We note that the set of {k| ^} are actually 
the discrete momenta for periodic boundary conditions. 
We will denote F m {k l ) = F m (kf ) + q) as F™(kf } ). It 
is easy to find that the many-body wavefunction of the 
ground state in the momentum space $ q ({kj }) is the 
Slater determinant of F^(k[ 0) ). By means of the sub- 
stitution gf- -> and -> the Chern 
number is derived to be 

C = -L / dq(V q $ q | x |V q $ q ) (3) 

2m JRq 

where the inner product includes a summation over 
{k<W}, and i? q denotes the rectangle of [0, f^] x [0, f^]. 

Using the Stokes theorem, we rewrite Eq. Q as a line 
integral 

C=^-i dl q ■ ($ q |V q $ q ) , (4) 

around the boundary of i? q , denoted by 9i? q . The Chern 
number given in Eq. Q is expressed as a winding number 
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FIG. 1: (Color online) The rectangle i? q has a size 4-7T 2 /L x L y , 
whose four vertices are denoted by q Q (a = 0, 1, 2, 3 ). 

along closed path dRq. We can divide dRq into N small 
line segments with q a (a = 0, 1 • • • N) as their endpoints. 
In each segment, one can replace the derivatives in Eq. (|4]) 
by discrete differences and the integral by a summation, 
yielding 

N-l 

C= — Arg[det(C Q , Q+1 )] , (5) 

where Arg(-) stands for the principal argument, and 
C a , a +i is a M x M coupling matrix^, with elements 

C^2+i = i F qJ F ^+i)- Here > we have used the relation 
(<frq a |$q Q+1 ) = det(C a . a+ i). Equivalently, one can first 

multiply the coupling matrices C = Yl^=o C a , a +i, and 
then diagonalize C. The Chern number is given by the 
sum of the phases of the eigenvalues of C divided by 27r. 

In practice, it is sufficient to take N — 4 and q Q with 
a = 0, 1, 2, and 3 (q4 = qo) to be the four vertices of i? q 
shown in Fig. [1] when L x ^> 1 and L y ^> 1. It is inter- 
esting to notice that for these q a , k'°) + q a still belong 
to the set of {k( )}, so all the quantities that are needed 
for the calculation of the Chern number can be obtained 
in the system with periodic boundary conditions. More- 
over, using the inverse FT, the matrix elements of the 
coupling matrix C a , a +i can be expressed in real space as 

CT a+ i = (^o|e l(qQ - q ° +l) - r |^ = o) • (6) 

Here, <p™ = q{v) with m — 0, 1, 2 • • • (M — 1) are the single- 
particle wave functions of the occupied electron states 
for periodic boundary conditions. The product of the 
coupling matrices 

C = Co,iCi, 202,3^,0 , (7) 

is first diagonalizcd to find M eigenvalues denoted as A m . 
Since C is not Hamiltian, its eigenvalues A m are usually 
complex numbers. Then the Chern number is given by 

1 M-l 

C = ^ E Ar §( A ") . ( 8 ) 

m=0 

where the range of the principal argument function Arg(-) 
is taken to be (— w, ir]. We mention that C tends to be 
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diagonal when M is sufficiently large^, so the sum of the 
phases of all the diagonal elements in C divided by 2-7T 
yields a good approximation to the Chern number. 

Prodan, Hughes, and Bernevig^i proposed an efficient 
numerical approach to the Chern number, by extending 
the formula for pure systems straightforwardly to disor- 
dered systems. Their expression is also in real space, 
which however involves multiple commutators of the co- 
ordinate operators and projection operator, and is un- 
likely reducible to a single term. Our method has the 
advantage that the formula contains only a single term 
and is more physically transparent. In the next two sec- 
tions, numerical calculations will be carried out to show 
the validity and efficiency of this formula for calculating 
the Chern number. 



III. THE HALDANE MODEL 

Let us consider the Haldane model 3 for the quantum 
anomalous Hall effect defined on a honeycomb lattice, 
including a random on-site disorder potential 

h = - d t d i +it v a4 £ j + Y w * 6 ^ ■ ( 9 ) 

Here, the first term describes the usual nearest-neighbor 
hopping with the hopping integral being taken to be 
unity, and the second term stands for the next-nearest- 
neighbor hopping with a complex hopping integral and 
Vij = (dkj x d lk ) z /\(d kj x d lk ) z \. i and j are two next 
nearest neighbor sites and k their common nearest neigh- 
bor. The vector di k points from k to i, and the distance 
between two nearest neighbor sites is taken to be unity. 
The third term represents a random on-site potential, 
with u}{ uniformly distributed between [W/2, W/2]. It 
is known 3 , that in the absence of disorder, the energy 
spectrum of the Haldane model has a middle band gap 
between — 3y/3t and 3y/3t. When the Fermi energy lies 
in the band gap, the system shows a quantized Hall con- 
ductivity e /h, without an applied magnetic field. 

Numerical calculation of the Chern number is carried 
out for a disordered system with size N=48 x 48 and 
t = 0.2, and the results are shown in Fig. [2] In the ab- 
sence of disorder (W — 0), the Chern number is well 
quantized to 1 within the middle band gap, which is con- 
sistent with the known result 3 -. With increasing disorder 
strength, the C = 1 plateau narrows, manifesting the 
levitation and pair annihilation^ of extended states for 
the valence and conduction bands, whose Chern numbers 
have opposite signs. The calculated Chern number at a 
fixed Fermi energy (E = 0) for t = 0.1 as a function of 
disorder strength is plotted in Fig. [3] for different sample 
sizes. It is found that the Chern number is robust against 
weak disorder W < 4. With increasing W from 4 to 6, 
the Chern number continuously decreases to nearly zero. 
With increasing the sample size, the transition process 
becomes sharper and sharper, which conforms the expec- 
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FIG. 2: (Color online) Chern number for the Haldane model 
as a function of electron Fermi energy E for t = 0.2, for several 
different disorder strengths. The result is averaged over 400 
disorder configurations for sample size N — 48 x 48. 
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FIG. 3: (Color online) Chern number for the Haldane model 
as a function of disorder strength W at t = 0.1, for several 
different sample sizes. The result is averaged over 400 disorder 
configurations. 



tation that the transition should become a sudden drop 
from 1 to in the thermodynamic limit. 



IV. HOFSTADTER MODEL 

We next turn to the Hofstadter models on a 2D square 
lattice. Under an external uniform magnetic field with 
the Landau gauge A — (yB, 0, 0), the model Hamiltonian 
is given by 

#o = - 53(e iw *4 >n <w fll „ + cl n c m , n+1 ) + H.c, (10) 

m , n 
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FIG. 4: (Color online) Calculated Chern number and elec- 
tron density of states for the lattice Hofstadter model in the 
full energy band for magnetic flux q/p — 1/16. The disorder 
strength is set to be W = and N — 64 x 64. The inset is a 
drawing of partial enlargement. 
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FIG. 5: (Color online) Chern number as a function of disorder 
strength W for the lattice Hofstadter for magnetic flux q/p = 
1/16 and sample size N = 64 x 64. The result is averaged 
over 400 disorder configurations. 



where integers m and n are the x and y coordinates of a 
lattice site, c m , n is the fermion annihilation operator on 
the site, and the hopping integral has been set to be the 
unit of energy. The magnetic flux quanta per plaquette is 
4> = 2nq/p. In the absence of disorder, the original band 
splits into p Landau levels for mutually prime integers 
p and q. The special case of q = 1 has been widely 
studied in previous works.— ~ 17 ' 26 ~ — It was found that 
each Landau band carries a Chern number +1 except 
the band at the center (E = 0) which carries a negative 
Chern number — {p — 1) for odd p or — (p/2 — 1) for even 
p. Figure 0] shows the calculated Chern number as a 
function of the Fermi energy at q/p = 1/16 for a sample 
of size N=64 x 64. We see that the Chern number of 
each Landau level is +1, except the central two Landau 
levels, each one having a Chern number —7. As a result, 
the sum of the Chern numbers of all the Landau levels 
in the full energy band is zero. 

To study the disorder effect, we include into the 
Hamiltonian a random on-site potential of the form 
E m ,„ w mi „cj„ in c min , with w m<n randomly distributed be- 
tween [— W/2,W/2]. The calculated Chern number at a 
fixed Fermi energy E = —2.75 is displayed as a function 
of disorder strength W in Fig. [5j The Chern number is 
well quantized for weak disorder. With increasing W , the 
interesting direct transition^ from the C = 2 quantum 
Hall plateau to a trivial insulator state with C = is 
reproduced. Moreover, the critical disorder strength is 
estimated to be about W — 3.5, in good agreement with 
the result obtained from the Kubo formula calculation^. 



V. SUMMARY 

To conclude, we have proposed an efficient coupling- 
matrix method for calculating the Chern number of dis- 



ordered systems. The present approach is applied to both 
the Haldane model and lattice Hofstadter model. The 
calculated Chern number is found well quantized even if 
the sample size is not very large. The calculated results, 
in particular the disorder-induced phase transition, are in 
good agreement with the known results. Our approach 
can be directly applied to calculate the spin Chern num- 
ber a n i 12 for the quantum spin Hall systems, from which 
the Z2 index can be extracted. 
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